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Abstract 

For one-qubit pure quantum states, it is already proved 
that the Voronoi diagrams with respect to two distances — 
Euclidean distance and the quantum divergence — coin- 
cide. This fact is a support for a known method to calculate 
the Holevo capacity. To consider an applicability of this 
method to quantum states of a higher level system, it is es- 
sential to check if the coincidence of the Voronoi diagrams 
also occurs. In this paper, we show a negative result for 
that expectation. In other words, we mathematically prove 
that those diagrams no longer coincide in a higher dimen- 
sion. That indicates that the method used in one-qubit case 
to calculate the Holevo capacity might not be effective in a 
higher dimension. 



1. Introduction 

The movement of trying to apply quantum mechanics 
to information processing has given vast research fields in 
computer science |6|. Especially among them quantum in- 
formation theory is developed as one of the richest research 
fields. Some aspect of quantum information theory is to in- 
vestigate a kind of distance between two different quantum 
states. Depending on the situation, several distances are de- 
fined in quantum states. 



The significance of quantum information theory is also 
due to its rich variety of mathematical methods used there. 
Especially, geometric interpretation of quantum states has 
been an important theme. The researches from that kind of 
view are categorized as "quantum information geometry." 
Actually some properties of a space of quantum states as a 
metric space have been researched in some contexts |3] |5J 
1. 

A quantum channel is a channel that transfers quantum 
information. Mathematically it is represented as an affine 
transformation between two Hilbert spaces, each of which 
is a representation of quantum states. One of the main prob- 
lem in the quantum information theory is how well a quan- 
tum channel transfers information. Especially, one of the 
most important indications that show a capacity of quantum 
channel is the Holevo capacity \ 2\. Intuitively speaking, the 
Holevo capacity indicates how much the channel preserves 
the size of the space of the quantum states. The measure 
used here is called the quantum divergence, which is a kind 
of a distance of two quantum states. 

Oto et al. |7| showed that the Holevo capacity of one- 
qubit states can be numerically computed by considering 
the image of the source points of the channel. In that paper, 
points in pure states are plotted so that they are almost uni- 
formly distributed with respect to the Euclidean distance. 
Here pure states are important because they appear in the 
boundary of the convex object which corresponds to the 
whole space of the quantum states. Although their images 



are dealt in a context of the divergence, the algorithm is 
reasonable because the spaces of pure quantum states with 
respect to the two distances have the same structure. In 
other words, in pure states, uniformly distributed points in 
the world of the Euclidean distance is also uniformly dis- 
tributed with respect to the divergence. The authors showed 
that fact mathematically considering Voronoi diagrams 1 3 1. 

A natural question that arises after this story is "What 
happens in a higher level systems?" If you could say the 
same thing in a higher dimension, the method used in the 
one-qubit case for a calculation of the Holevo capacity 
might be applied to a general case. However, unfortunately 
we found it is not the case. The main result of this paper is 
the fact that in a higher dimension, the two distance spaces 
— the space of pure states with respect to the Euclidean 
distance and the one with respect to the divergence — no 
longer have the same structure. More correctly, we prove 
that the Voronoi diagrams with respect to the two distances 
do not coincide in a higher dimension. Additionally we give 
some examples for the understanding of the structure. 

The rest of this paper is organized as follows. First in 
Section[2] we give some basic facts in quantum information 
theory. In Section[5] we explain briefly the known fact about 
one-qubit quantum space. Sectionals the main part of this 
paper, where we show some workout for a higher level case 
and give some illustrative examples. Lastly in Section|5l we 
summarize the result and give it some discussion. 

2. Preliminaries 
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Actually, any matrix which is Hermitian and whose trace 
is one is expressed this way with some adequate {£«}. 
This condition doesn't contain a consideration for a semi- 
positivity. To add the condition for a semi-positivity, it is 
not simple as in one-qubit case, and we have to consider 
complicated inequalities QJQ|. Note that this is not the only 
way to parameterize all the density matrices, but it is reason- 
ably natural way because it is natural extension of one-qubit 
case and has a special symmetry. 

Additionally our interest is a pure state. A pure state is 
expressed by a density matrix whose rank is one. A density 
matrix which is not pure is called a mixed state. A pure 
state has a special meaning in quantum information theory 
and also has a geometrically special meaning because it is 
on the boundary of the convex object. In one-qubit case, the 
condition for p to be pure is 



2.1. Parameterization of quantum states 

In quantum information theory, a density matrix is repre- 
sentation of some probabilistic distribution of states of par- 
ticles. A density matrix is expressed as a complex matrix 
which satisfies three conditions: a) It is Hermitian, b) the 
trace of it is one, and c) it must be semi-positive definite. 
We denote by S(C d ) the space of all density matrices of 
size d x d. It is called "d-level system." 

Especially in two-level system, which is often called 
"one-qubit system", the conditions above are equivalently 
expressed as 
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The parameterized matrix correspond to the conditions a) 
and b), and the inequality correspond to the condition c). 

There have been some attempt to extend this Bloch ball 
expression to a higher level system. A matrix which satis- 
fies only first two condition, Hermitianness and unity of its 
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This is a surface of a Bloch ball. On the other hand, in 
general case, the condition for pureness is again expressed 
by complicated inequalities. 

2.2. The quantum divergence and the 
Holevo capacity 

We define the log of density matrix. When eigenvalues 
of p are diagonalized with a unitary matrix X as 



p = X 
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the log of p is defined as 
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The quantum divergence is one of measures that show 
the difference of two quantum states. The quantum diver- 
gence of the two states a and p is defined as 



D(a\\p) =Trcr(loga-logp). 



(6) 



Note that though this has some distance-like properties, it 
is not commutative, i.e. £)(<t||p) ^ D(p\\a). The diver- 
gence D(a\\p) is not defined when p does not has a full 
rank, while a can be non-full rank. This is because for a 
non-full rank matrix, a log of zero appears in the definition 
of the divergence. However, since log is naturally de- 
fined as 0, some eigenvalues of a can be zero. 

A quantum channel is the linear transform that maps 
quantum states to quantum states. In other words, a linear 
transform T : M(C; d) — > M(C; d) is a quantum channel if 
Y(S(C d )) C S(C d ). 

The Holevo capacity |2| of this quantum channel is 
known to be equal to the maximum divergence from the 
center to a given point and the radius of the smallest enclos- 
ing ball. The Holevo capacity C(T) of a 1-qubit quantum 
channel T is defined as 

C(r)= inf sup D(T(a)\\T(p)). (7) 

3. One-qubit case 



necessarily uniform with respect to the divergence. So why 
does this mechanism work correctly? You cannot say the 
approximation is good enough unless uniformness of the 
image of points is guaranteed. That concern is overcome by 
comparing the Voronoi diagrams. The following theorem is 
more precise description 1 3 1. 

Theorem 1. Suppose that n one-qubit pure states are given, 
the following Voronoi diagrams of them are equivalent: 

1. The Voronoi diagram in pure states obtained by taking 
a limit of the diagram with respect to the divergence 

2. The Voronoi diagram on the sphere with respect to the 
ordinary geodetic distance 

3. The section of the three-dimensional Euclidean 
Voronoi diagram with the sphere 

Note that although the divergence is not defined in the 
pure states, we can consider the limit of the diagram with 
respect to the divergence in the pure states. In this paper, 
Voronoi sites are plotted as a first argument of the diver- 
gence D(-||-). From now on, when we say "Voronoi dia- 
gram with respect to the divergence", it means a limit of a 
diagram with sites in the first argument of 



Our first motivation to investigate a Voronoi diagram in 
quantum states is the numerical calculation of the Holevo 
capacity for one-qubit quantum states Q. In that paper, 
some points are plotted in the source of channel, and it is 
assumed that just thinking of the images of plotted points is 
enough for approximation. Actually, the Holevo capacity is 
reasonably approximated taking the smallest enclosing ball 
of the images of the points. More precisely, the procedure 
for the approximation is the following: 

1 . Plot equally distributed points on the Bloch ball which 
is the source of the channel in problem. 

2. Map all the plotted points by the channel. 

3. Compute the smallest enclosing ball of the image with 
respect to the divergence. Its radius is the Holevo ca- 
pacity. 

In this procedure, step 3 uses a farthest Voronoi diagram. 
That is the essential part to make this algorithm effective 
because Voronoi diagram is the known fastest tool to seek a 
center of a smallest enclosing ball of points. 

However, when you think about the effectiveness of this 
algorithm, there might arise a question about its reason- 
ableness. Since the Euclidean distance and the divergence 
are completely different, Euclideanly uniform points are not 



4. Higher level case 

In this section, we show that the coincidence which hap- 
pens in one-qubit case never occurs in a higher level case. 
To show it, it is enough to look at some section of the dia- 
grams with some hyperplain. If the diagrams do not coin- 
cide in the section, you can say they are different. 

Suppose that d > 3 and that the space of general quan- 
tum states is expressed as Equation (0, and let us think the 
section of it with a hyperplain: 



£,d+2 — £,d+3 — ■ ■ 

Then the section is expressed as: 
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The elements of this matrix are except diagonal, (0, 1), and 
(1,0) elements. This matrix is diagonalized with a unitary 



matrix as: 
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Now we will figure out the necessary and sufficient con- 
dition for the diagonal matrix of Equation d 1 Oft to be rank 1 . 
For that condition to hold, the following three cases can be 
considered: 



Case 1 (only e?-th raw of the matrix is non-zero) 

£1 = £2 = • • • = £d-i = -l, £d = £d+i = 0. 

Case 2 (only one i-th raw {3 < i < d — 1) is non-zero) 

£1=6 = -Urf = £ rf +i =0, 

all of £j (3 < j < d — 1) are —1 except one (let its index to 
be k) and £/. = d — 3. 



£i + £ 2 =d-2,^i + ^(d + d +1 ) = i, 

£3 = £ 4 = • • • = £d-i = -l. (17) 

Note that it is impossible that only A x is non-zero. In both 
Case 1 and Case 2, the set of points that satisfies the condi- 
tion is just one point, so our main interest is Case 3. The set 
of points that satisfies this condition is a manifold. Actually, 
Case 3 satisfies 

^^4<g + S + ,) = l, C>8> 

and this is an ellipsoid. 

Then we prepare for workout of the divergence. The log 
of p is expressed as: 
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Thus, we obtain 
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With some workout, we get 
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Using these fact and the assumption 771 + 772 = £1 + £ 2 
rf — 2, we get 



Tr a log p 



2r "I d 3 ? 



ilogAiA 2 . 



(22) 



Next we think of a Voronoi diagram with only two re- 
gions for simplicity. It is enough for our objective. Let 
a and a be two sites, and suppose that p moves along the 
boundary of the Voronoi regions. Suppose that a and a are 
parameterized by {rjj} and {fjj} respectively in the same 
way as p. 

We consider what happens if r(0 < r < 1) is fixed and 
the following holds: 



ft + 6 = d - 2, & 
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The condition < r < 1 means that p is semi-positive 
and not a pure state while r = 1 in pure states. In other 
words, we regard that p is on the same ellipsoid obtained by 
shrinking the ellipsoid expressed by Equation ( 1181 . These 
settings are in order to take a limit of a diagram to get a 
diagram in the pure states. Taking the limit r — > 1, we can 
get a condition for pure states. This procedure is analogous 
to the method used in 1 3 1 . 

Now to think of the shape of boundary, we have to solve 
the equation 



D(a\\p) = D(a\\p), 
and this is equivalent to 

Tr (o--cr)logp = 0. 
Using Equation J22i . we obtain 

Tr (tr — cr) log p = 

(Vd-ijd)£d + (%+1-%+iKd+i 



(24) 



(25) 



2r 



4(ru^!)(6-^) 



d 2 



(26) 



Here when r = 0, this is zero because A1/A2 = 1. In 
that case, p can take only one point, but we do not have to 
care about this case because we are going to take the limit 
r — > 1. From now on, we suppose r > and that means 
Ai/A 2 ? 1. 

Hence we get the following equation that holds in the 
boundary of the Voronoi diagram: 
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Consequently, taking the limit r — > 1, we get Equation ( I27i 
as the expression of the boundary in pure states. 

A careful inspection of Equation (127 \ tells us a geomet- 
ric interpretation of this boundary. We obtain the following 
theorem: 

Theorem 2. On the ellipsoid of the pure states which ap- 
pears in the section with the hyperplain defined above, if 



transfered by a linear transform which maps the ellipsoid 
to a sphere, the Voronoi diagram with respect to the diver- 
gence coincides with the one with respect to the geodesic 
distance. 

Proof. Think of the affine transform defined by 



*\ 




/ c d ~ 2 
1 £1 — 

d 


y 


- 


id 


* ) 




\ Cd+l 



(28) 

then Equation (127 \ is expressed as 

x (x — x) + y (y — y) + z (z — z) = 0, (29) 

while Equation Jl 81 becomes 

x 2 +y 2 + z 2 = l. (30) 

Thus when (x, y, z) and (x,y,z) are fixed, the point 
(x' , y', z') which stand for 77 runs along the geodesic. □ 

Now we work out the Voronoi diagram with respect to 
Euclidean distance. Under the assumption above, the Eu- 
clidean distance is expressed as 

d(a, p) 

^{vi-iifHrn-^f+im-^fHrn+i-^d+i? 

= 2(77! - a) 2 + fad - id) 2 + (Vd+i - td+i) 2 , (31) 
and we get the equation for boundary as 

d(a,p) - d(a,p) = 

- 4(771 - 771)^1 - 2(ij d - fjd)£d - 2(77^+1 - fjd+i)£d+i 
+ 2(77? - it) + ( v 2 - f, 2 ) + ( V 2 d+1 - ff d+1 ) = 0. (32) 

By comparing the coefficients of £1, £<j> and £d+i, we can 
tell that the boundaries expressed by Equation i21l and i32\ 
are different. To show how different they are, we give some 
examples in the rest of this section. 

Example 1. Suppose that (r]i,r)d, Vd+i) = (d — 1, 0, 0) and 
{VuVd,fjd+i) = (-1,0,0), then the boundary is a) £1 = 
for the divergence, and b) £1 = 1 for the Euclidean 
distance. Fig. [2 shows this example for d — 5. In Fig. [2 
Voronoi sites are located on the top and the bottom of the 
ellipsoid. The two diagrams are the same when d = 4, but 
are different otherwise. 

Example 2. Suppose that (771,77^,77^+1) = (0,1,0) and 
(771,77^,77^+1) = (0,-1,0), then the boundary is, for 
both the divergence and Euclidean distance, expressed by 

6i+i = o. 




Figure 1. An example of a Voronoi diagram 
with two sites. The figure on the left is the 
diagram by the divergence, and the figure on 
the right is the diagram by the Euclidean dis- 
tance. 



Example 3. Consider the Voronoi diagram with the follow- 
ing eight sites: 

f d-2 d J_ J_\ 
\ 2 + 2V3' V3' 

fd-2 d [2 \ 

fd-2 d [2 \ 

where ±'s mean all the possible combinations. Then the 
Voronoi diagrams look like Fig. [2] This figure is also for 
d = 5. Obviously they are different. 




Figure 2. An example of a Voronoi diagram 
with eight sites. The left is the diagram by the 
divergence, and the right is by the Euclidean 
distance. 



expression of some section of the boundary of the Voronoi 
diagram with respect to the divergence. The section is an 
ellipsoid and after some linear transform, the boundary be- 
comes a geodesic on a sphere. Interestingly this is similar to 
the whole space of the one-qubit states even though a space 
of higher level has much more complicated geometric struc- 
ture. 

We also showed some geometric structure concerning 
how pure states appear in a whole quantum states. Although 
our result is very restricted, we believe this will be a help for 
further understanding of the structure. To relax the restric- 
tion is our future work. 

The result shown in this paper depends on the parame- 
terization of density matrix. The parameterization used in 
this paper, though it is very natural one, is not unique. To 
think of another parameterization is another future work. 
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5. Conclusion 

We proved that in n-level system for n > 3, the Voronoi 
diagrams with respect to the divergence and Euclidean dis- 
tance do not coincide. Additionally we obtained an explicit 



